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Abstract 
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Introduction 

Recent neutron star observations make it possible to distinguish 
in the entire population of neutron stars a class of massive stars. 
Particularly important data that indicate the existence of massive 
neutron stars concern the 3.1 ms radio pulsar in the binary system 
J1614-2230 [1] and the binary pulsar J0348-H0432 [2] with the 
rotation period P = 39. 1 ms. These binary systems possess features 
that make them extremely important for astrophysics. PSR J0348-I- 
0432 includes a neutron star and a low mass white dwarf 
companion (M jy£) = 0.172MQ), the orbital period of this system is 
2.4 hr. The estimation of the neutron star mass based on radio 
timing and the precise spectroscopy of the white dwarf companion 
gives as a result Mi^s = 2.01 + 0.04Mq . This binary system due to 
its special characteristics, namely: the neutron star and the low- 
mass white dwarf and tight, relativistic orbit is expected to provide 
information on the orbital decay due to the emission of 
gravitational waves. PSR J1614-2230 is a binary system that 
consists of a neutron star and a white dwarf. Measurements of the 
general relativistic time delay (Shapiro effect) allows for very 
accurate estimation of the mass of both components: neutron star 
mass Af/v^s = 1-97 + 0.0471^0 and white dwarf mass 
MwD = 0.500 + 0.006Mq . 

In the light of these observations, it is important to know the 
evolutionary path that leads to the formation of a neutron star with 
such large mass. The most general scenario predicts that the 
formation of a cold catalysed object is preceded by a very specific 
phase of a hot neutron star with trapped neutrinos, which defines a 
proto-neutron star [3,4]. Understanding the evolution of neutron 
stars requires a detailed analysis of the properties of a proto- 
neutron star matter, especially the heat transport [5]. 

The existence of massive neutron stars entails important 
consequences for the equation of state (EoS) of dense nuclear 
matter. In general, the precise measurements of pulsar masses put 
constraints on the form of the EoS and make problematic the 



existence of exotic particles such as hyperons and quarks [6] in the 
very inner part of a neutron star. Realistic models of neutron stars 
provide stratification of their internal structure by distinguishing 
two main components: a crust and a core. The crust, that splits 
into the outer and inner part, describes the outer layer of a neutron 
star with subsaturation densities and contains only a small 
percentage of a neutron star mass [7]. The description of the 
core of a neutron star is modelled on the basis of the EoS of dense 
nuclear matter in a neutron-rich environment [8] having a density 
that ranges from a few times the saturation density (mq) to about an 
order of a magnitude higher and at such densities hyperons are 
expected to emerge [9]. 

Despite the fact that neutron star matter is directly affected by 
the nature of strong interactions, it is not possible to give its 
description on the basis of quantum chromodynamics (QCD) even 
though it is fundamental theory of strong interactions. The 
description of nuclear matter is based on different models and with 
the use of different approaches. The commonly adopted method is 
the mean field approximation either relativistic or nonrelativistic. 
In the latter case the Skyrme forces [10], which contain 
parameters that have been established by adjusting nuclear matter 
and finite nuclei properties, yields very good results. In the 
microscopic approach the models of the nuclear matter EoS has 
been constructed with the use of the Brueckner-Hartree-Fock 
(BHF) approximation with the employed realistic two-nucleon 
interactions. In order to obtain the correct description of nuclear 
matter properties a phenomenological three-nucleon interactions 
was introduced [11,12]. The paper [13] comprises a first attempt 
to describe a properties of strangeness-rich matter within 
Brueckner theory. The extension of the analysis for the ;8-stable 
nuclear matter with hyperons was given in the paper [14]. Further 
development of the method that allows one to model strangeness- 
rich nuclear matter in this microscopic approach takes into 
account not only hyperon-nucleon (YN) but also hyperon-hyperon 
(YY) interactions [15]. In the context of massive neutron stars that 
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contain hyperons the efTect of the hyperonic three-body forces was 
estimated [16]. Calculations of the EoS performed on the basis of 
BHF approximation with additional phenomenological density- 
dependent contact terms do not result in a correspondingly large 
neutron star masses but form the basis for further more 
sophisticated analysis. 

Observations of neutron stars with masses of the order of 2Mq 
also do not excluded the existence of hybrid neutron stars with 
quarks nucleated in their inner cores. Elforts were done to analyse 
conditions under which the appearance of a quark core does not 
cause too strong softening of the EoS. There were several attempts 
to construct models of hybrid stars that satisfactorily reproduce the 
mass of a neutron star at the level of 2Mq. The differences 
between them result from the form of the EoSs that were used to 
describe the hadronic and the quark phases. Models that were 
used to describe the hadronic phase are the phenomenological 
models with the most popular relativistic mean field models [17] 
with cUfiFerent parameterizations and the microscopic Brueckner- 
Hartree-Fock approach [18-20]. The deconfined quark phase was 
modelled for example on the basis of the MIT bag model or the 
Nambu-Jona-Lasinio (NJL) model [20-22]. However, it was found 
that the transition to the quark matter phase is possible only for 
selected EoSs and parameterizations. Detailed analysis performed 
for the model of a hybrid star constructed on the basis of 
microscopic Brueckner-Hartree-Fock approach for hadronic 
phase and the MIT bag model, the IS^L model and chromodi- 
electric model for the quark phase are given in the following 
papers [18-20]. Also the effect of a hyperonic three-body force on 
the metastibiUty of compact stars was investigated leading to the 
estimations of the maximum mass of a hybrid star at the level of 
1.62Mq [20]. Results that were obtained for the hybrid star 
model with the quark matter phase described in the framework of 
a standard color superconducting NJL model and the hadronic 
phase constructed on the basis of the Dirac-Brueckner-Hartree- 
Fock EoS for the Bonn-A potential indicate the possible existence 
of massive hybrid stars [2.'5]. 

The purpose of this paper is to study the impact of hyperons on 
the properties of the matter of neutron stars, and on their 
structure. In general, analysis of the role of strangeness in nuclear 
structure in the aspect of multi-strange system is of great 
importance for both nuclear physics and for astrophysics, and 
leads to a proper understanding of the properties of a hyperon star. 
The aim of this paper is an accurate analysis of the properties of 
neutron stars with particular regard to their internal structure. The 
key issue to be examined relates the determination of the factor 
which in the considered model is responsible for the stiffening of 
the EoS. Another aim of this work is to provide understanding on 
the impact of the model with the extended vector meson sector on 
the internal structure of a neutron star. 

The relativistic approach to the description of nuclear matter 
developed by Walecka [24] is very successful in describing a 
variety of the ground state properties of finite nuclei. Although the 
original Wale[:ka model properly describes the saturation point 
and the data for finite nuclei, it has been insufficient to properly 
describe the compression modulus of symmetric nuclear matter at 
saturation density. The nonfinear self-interactions of the scalar 
field (the cubic and quartic terms) were added in order to get an 
acceptable value of the compression modulus [2.5,26]. Additionally 
the inclusion of a quartic vector self-interaction term softens the 
high density component of the EoS [27]. The estimation of the 
incompressibflity coefficient of symmetric nuclear matter Kq, 
which is made on the basis of recent experimental data, points to 
the range 240+ 10 MeV [28]. 



Models that satisfactorily reproduce the saturation properties of 
symmetric nuclear matter lead to considerable differences in a case 
in which asymmetry dependence is included [29,30]. Thus, the 
proper model of the matter of neutron stars requires taking the 
effect of neutron-proton asymmetry into consideration. This, in 
turn, leads to the inclusion of the isovector meson p. The standard 
version of the introduction of the p meson field is of a minimal 
type without any nonlinearities. This case has been further 
enlarged by the nonlinear mixed isoscalar-isovector couplings, 
which modify the density dependence of the p mean field and the 
symmetry energy [31-33]. The analysis of the nonlinear models 
should include results obtained for the relativistic FSUGold 
parametrisation [34]. 

A theoretical description of strangeness-rich nuclear matter 
requires the extension of the model to the full octet of baryons and 
additional meson fields were introduced to reproduce the 
hyperon-hyperon interaction. The model that is considered is 
constructed on the basis of the hadronic SU{3) theory, which 
naturally includes nonlinear scalar and vector interaction terms. 
The characteristic feature of the model is the very special form of 
the vector meson sector, which permits more accurate description 
of asymmetric strangeness-rich neutron star matter [35]. The 
primary goal of this paper is to maximize the understanding of the 
influence of the nonlinear vector meson couplings on the form of 
the EoS and through this on a neutron star structure. 

Having obtained the EoS the analysis of the maximum 
achievable neutron star mass for a given class of models can be 
performed. Observational results limit the value of a neutron star 
mass and thereby put constraints on the EoS of high density 
nuclear matter. Recent observations point to the existence of a 
high maximum neutron star mass [1,2], what is inconsistent with 
theoretical models that involve hyperons. 

The model 

Recent observations of the binary millisecond pulsars J 1 6 1 4- 
2230 [1] and J0348+0432 [2] have led to the precise estimation of 
neutron star masses: (1.97 + 0.04)Mo and (2.01 ±0.04)Mo. This 
places the maximum neutron star mass at rather high values and 
rules out most of the EoSs with hyperons as models that involve 
exotic particles predict maximum neutron star masses well below 
the stated values. There is a need to analyse whether it is possible 
to construct an EoS of neutron star matter that gives adequately 
high maximum mass despite including hyperons. 

In this paper a description of nuclear matter based on an 
effective model constructed within the framework of the nonlinear 
realization of the chiral SU(3)i^x SU(i)ji symmetry [36-38] is 
giving. Details can be found in the papers [35,39]. Baryons and 
mesons constitute the basic degrees of freedom of the model, and 
consequentiy the Lagrange density function £ splits into parts that 
are adequate to describe baryon Cb and meson Cm sectors 
supplemented by the term that represents baryon-meson interac- 
tions £,„;, and takes the form C = Cg-'r Cm+ C-im- The meson 
sector of the considered model includes spin zero and spin one 
meson states. Nonets of cfifferent meson types, spin zero (scalar) 
and spin one (vector), can be written as the sum of the singlet and 
octet matrixes M =Msm -\-Moct- Under the assumption of SU(3) 
symmetry, a very general form of the interaction Lagrangian 
includes a mixture of the symmetric [D -type) and antysymmetric 
{F -type) couplings and the S -type coupling that denotes the 
meson singlet state interaction 
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+ (1 - ^m)\BBM\^) - ^ [BBMU 

= - Vlg^ [ocM(Tr(BMoc,B) ~ Tr{BBMoc>)) 
+ (1 - aM){Tr{BMoc,B) + Tr(BBMoc,))\ 
1 

7^' 



-gf^Tr(BB)Tr{Msi„), 



(1) 



where the expUcitly given baryon matrix B has the form 



The <l> meson taken as pure ss state leads to the ideal mixing 
with 0]7«35.3°. If the determination of the baryon-vector meson 
couplings bases on the assumption that nucleons do not couple 
to the ss meson, then g]' = y^gg- In the limit ay = l (only the 
F -type couphng remains), the couphng constants are related to 
the additive quark model and the vector meson coupling constants 
are given by the following relations: 



gAco = g-Lio = 2gSco = 2 SNco = 2g8 , 



gA^=g-£4-- 



gs4 _\/2 



gNo>- 
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The high density limit of the EoS of dense matter in neutron 
star interiors is dominated by the contribution that come from the 
baryon number density, which justifies the construction of a model 
(2) that includes a broad spectrum of mixed vector meson coupKngs. 
The extended vector meson sector, which stems from the SU(3) 
invariants, can be written in the following form 



Generally, the baryon-meson interaction is characterised by the 
following coupling constants: the octet g^ and singlet gj^ 
couphng constant, the parameter aM, which stands for the 
F/(F + D) ratio, and a mixing angle (?,v( that relates the physical 
meson fields to the pure octet and singlet states. The index M 
concerns the scalar (A1 = S) and vector {M = V) mesons 
(pseudoscalar and axial vector mesons, which have a vanishing 
expectation value at the mean field level, are not considered in this 
model). In the case of the scalar meson sector baryon masses are 
generated by the vacuum expectation value that is attained by two 
scalar meson fields and the parameters gi,gi, and as have been 
chosen to fit the experimental values of the baryon-octet masses 
[36,37]. Considering the vector meson sector of this model the 
octet matrix 



(3) 













Voc,= 


P~ 
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supplemented with the singlet state 
1 



Vsin = -^diag{yo,yo,vo) 



Cv=\c{Tr{VV)f+]^dTr{{VVf)+y{Tr{V))\ 



where V represents the matrix of the vector meson fields, and 
c,d/ are the general coefficients that have been determined by 
assuming that in the case of a neutron star with zero strangeness, 
the model described by the TMl parameter set [27] is recovered; 

C='^Ci-Av(gNpgNaf', d=-Ci-\-2Kv{gNpgNo,f, / = 0. This 

assumption leads to the Lagrangian function, which embodies 
contributions from the baryon and meson sectors supplemented by 
the parts that describe the baryon and meson interactions 



B ^ 

1 2 2 1 3 1 4 
-^^o^^ -gftf^ -4^4ff 



(6) 



comprises the vector meson nonet. The combinations of the 
unphysical SU(3) singlet (vq) and octet (vg) states produce the 
physical CO and (j) mesons 



(p = vs cos 6v — vo sin Oy, 
(0 = vg sin 9y + vq cos 6y. 



(4) 



where the covariant derivative equals Dn = 8fi + igBiiiCOn + 
igB<li<t'ii + ^SBp^BPfi, denotes isospin of baryon B. The baryon 
effective mass is defined as follows meff,B = fnB — gBa^ — gBrT'('*, 
while Q^v, R/iv, and <!),„, are the field tensors of the co, p, and (j) 
mesons. A proper description of hyperon— hyperon interaction 
requires the presence of hidden strangeness mesons: scalar (a*) and 
vector {(j)). The Lagrangian fimction (6) describes the )S-equilibrat- 
ed neutron star matter; thus there is also a need to consider the 
Lagrangian of free leptons £/ 
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(7) 



(10) 



All nonlinear vector meson couplings that occur in this model 
have been brought together in the form of a vector potential 



'*'0-2c3 Wo =S'5«^G«*' (11) 



1 



1 



)=-C3(co"a;,)^ + -C3(p''V;)' 

+ AF(g^a,g^,)'(«''«,)(p''VP 
1 /I 



+ )^{^Ct,- kvigNojgNp)' 



(8) 



n]ff^/o-2c^il=gBpdzni„ (12) 



'"<#,(/o - (c3 + 2AFfe«^Bp) )/o = gB^^snb, ( 1 3) 



(i-fd^-meffji-gBioii -XB)y'^wo 
-gBphBy'^-r'ro -gBit,XBV%)fl/B = 0. 



(14) 



The description of dense, hyperon-rich nuclear matter given by 
the Lagrangian (6) in the case of non-strange matter is reduced to 
the standard TMl model with an extended isovector sector. This 
extension refers to the presence of the m — p meson coupling and 
enables modification of the high densit)' limit of the symmetry 
energy. The strength of this coupling is characterised by parameter 
Ay. For each value of parameter Ay, the parameter gffp has to be 
adjusted to reproduce the symmetry energy E^y^ = 25.68 MeV at 
A:f = 1.15 fm-i [40]. 



where Jb and Ijb denote the spin and isospin projection of baryon 
B, meff,i, {i = co,p,<j>) are elfective masses assigned to vector meson 
fields: 



ff^lff, w='ni + 3C3wl+2Ay(gBo,gBpfrl 



+ ( ^ C3 - ArigBmgBpf )fo 



(15) 



The equation of state 

The mean field approach has been adopted to calculate the 
EoS. In this approximation, meson fields are separated into 
classical mean field values: so,SQ,wo,ro/o and quantum fluctua- 
tions, which are neglected in the ground state: 



'nlff.p =ml + 3c3rg -I- 2Ay(gBa>gBpfwl 



+ ( ^C3 -Ay(gBwgBpf )fo 



(16) 



(7 = ff + .so, a*=d*+sl, 

COf,=Wp + Wo, Wo =<COf,> Sop = < (Bo > , 

■;=j6;+'-o, ro=<p;;><5o„()'"=<Po>, 

<t'p=(l>ii+fo> h=«t>p>5op=<h>- 



The preferable attribute of the considered model is its very 
diverse vector meson sector, which allows one to study the 
relevance of different vector meson couplings for the form of the 
EoS. The Lagrangian function (6) makes it possible to calculate the 
equations of motion from the corresponding Euler-Lagrange 
equations. The obtained results, written in the mean field 
approximation, take the form: 



'lff,^=rr^^+[jC^-Ay(gB^gBpfywl + r^^ 



^Cj,+Ay{gBa>gBpf )/o^- 



(17) 



Analysis of the equations of motion raises the issue of the 
medium effects on the properties of the hadronic matter namely 
the effective baryon and vector meson masses. Considering the 
case in which nucleon mass depends only on non-strange 
condensate and referring to the Walecka model, the relation for 
the reduced effective baryon masses in the medium in the 
relativistic mean field description [31] can be formulated as 



meff,B = meff^{S(),sl)=mB-gBoSO-gBn*sl, 



(18) 



mlso=-gT.sl-gAsl 



, ^ 2Jb+\ fp'B 



mef)\B{so,sl)k^ 



where the terms gBa-^i) and gBa'S}) reprcsc'nt tlu' modification of 
(9) baryon masses due to the medium and niB denotes baryon masses 
generated by the vac:uum expectation values attained by scalar 
meson fields [36,37]. 

The source terms in equations (11-13) can be expressed by 
introducing the parameters &q and &s, which measure the 
contributions of strange {xb) and non-strange quarks: 
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Table 1. Vector meson coupling constants. 





Baryon (B) 


xb 


gB,„ 


gB^=XBg.\,:, 


gBp 


n 


0 




0 


gp 


P 


0 


g,>, 


0 


gp 


A 


1 

3 


2 


v/2 


0 


Z 


1 

3 


2 


V2 


2gp 




2 
3 


1 

3^'" 


-2^g,. 


gp 



Baryon-vector meson coupling constants, = (1 —-XB)gNoj and xb counts the contribution of strange quarks, gNoj^gw and g^^p =gp. 
doi:10.1371/journal.pone.0106368.t001 



Sq=J2(1-xb)-^, (19) 



S3=ThB-, (20) 

5s=y^x,'^. (21) 

Thus, the number density of u and d quarks refers to <5g, 
whereas the number density of strange quarks is given in terms of 
the parameter 5s- The parameter St, aims to evaluate the 
dilference between the density of u and d quarks, ni,= "^^riB 
denotes the total baryon number density. The baryon-vector 
meson coupling constants determined from the symmetry relations 
have been summarised in Table 1. The parameters 3q, 5s and 5t, 
wiU be further used to define the asymmetry parameter fa and 
strangeness content fs of the matter of a neutron star: 



where 5*^ denotes the strangeness of baryon B (5 = A,S,H). The 
numerical solution of the equations of motion, which depends on 
the form of the elfective vector potential t/jJoJJ/, has been limited to 
the result with only one single real solution. The demand for the 
existence of one real solution puts constraints on the value of 
parameter Ay. Such an analysis in the simpler case of symmetric 
nuclear matter leads to an equation that relates Ay and C3 
parameters and enables the critical value of parameter Ay^cr to be 
estimated 

- lcl+4Ayc,(g^g„f-AUg^g„f = 0, (24) 

where g„j=gNio and gp =gNp- 

The existence of one single real solution is not satisfied for 
Ay>Ay^cr- The Critical value of the parameter Ay^^r calculated 
for the selected parameterisations are collected in Table 2. In 
order to calculate the energy density and pressure of the nuclear 
matter, the energy momentum tensor T^,, which is given by the 
relation 



A=^. (22) ^ 8C 



rib 



^ has to be used. In equation (25) (9, denotes the boson and fermion 

fs==-^ , (23) fields. The energy density £ is equal to <7oo>, whereas the 

''^ pressure V is related to the statistical average of the trace of the 

spatial component Ty of the energy momentum tensor. Calcula- 



Table 2. The limiting values of the parameter Ay. 





Parameter set 


C3 


^v.c,■ 


K„ (MeV) 


NL3 [61] 


0 




271 


FSUGold [34] 


418.39 


0.0517 


230 


TMA [62] 


151.59 


0.0318 


318 


TMl* [63] 


134.624 


0.0215 


281.1 


TMl [27] 


71.3 


0.0156 


281.1 


TM2 [64] 


84.5318 


0.0186 


343.8 



The limiting value of the parameters Ai- ^, and cs together with the incompressibility of the symmetric nuclear matter Ka taken at the saturation density calculated for 

the chosen parameter sets. 

doi:l 0.1 371/journal.pone.Ol 06368.t002 
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tions done for the considered model lead to the foUowing explicit 
formulas for the energy density and pressure: 



+ \{^<:3+^vig<ogpf^Jt 

B "■ JO 



(26) 



This result agrees with theoretical predictions for H in nuclear 
matter obtained in the model D of the Nijmegen group [44]. 
However, the missing-mass spectra of a double-charge exchange 
reaction (K^ ,K^) on a '^C target have suggested the H well 
depth of 14-16 MeV (-16 MeV< t/^'^' < - 14 MeV) [45,46]. 
Data related to the 2 -nucleon interaction are controversial. 
Analysis of the experimental data indicate a repulsive S -nucleus 
potential [47], with a substantial isospin dependence [48]. In the 
case of YY interactions, the only sources of information are the 
double-strange hypernuclear systems. Several events have been 
identified that suggest an attractive AA interaction. The most 
promising results, known as the NAGARA event [49] with the 
^^He hypernucleus, indicate that the AA interaction is weakly 
attractive. The determined value of this potential at the level of 
= — 5 MeV permits a parameter set which reproduces this 
weaker AA interaction to be estimated [39]. 

The potential that describe hyperon-nucleon and h^'peron- 
hyperon interaction can be written in the form that involves both 
the scalar and vector coupling constants 



'P= \miyo+\mlrl+]^mlf^-]^mlsl 

^ 3 ^ 4 ^^*2,^/4,4\ 

+ ^rigrngpfw^rl + ^ C3 - Arigo^gpY) (wg + r^Yo 



(27) 



+ \(\c3+Av(go,gpf )/o' 



1 

3^ 



'F.B 



0dk 



0 



+in; 



where £l and Vl denote the contributions coming from leptons. 

Coupling constants 

Understanding the nature of interactions between baryons is a 
decisive factor for the properties of neutron stars especially in the 
case when hyperons are included. A precise description of baryon 
interactions in the strange sector of the model is essential for the 
correct construction of the EoS. However, the incompleteness of 
the experimental data intensifies the uncertainties that are 
connected with the evaluation of coupling constants that involve 
strange baryons. In general, there are very few data available to 
describe hyperon-nucleon (YN) and hyperon-hyperon (YY) 
interactions. 

Hyperon-vector meson coupling constants are taken from the 
quark model. They are summarised in Table 1. In the scalar 
sector, the scalar couplings g^a of the A, Z and E hyj)C'roiis rerjuire 
constraining in order to reproduce the estimated values of the 
potentials felt by a single A, 2 and H in the saturated nuclear 
matter. In the case of A hypernuclei there is a considerable 
amount of data on binding energies and single particle levels 
allowing the identification of the potential felt by a single A in 

nuclear matter in the range -30 MeV < t/^^' < - 27 MeV 
[41,42]. Considering the H hypernuclei current knowledge about 
the interaction of H hyperons with nuclei is very limited. Dover 
and Gal [43] based on early emulsion data indicated an attractive 
H -nucleus potential in the range —24 MeV <c4^<-21MeV. 



If 



(B) 



--gYaSQ -gYuiWQ +gYa*sl -gYijifo = 
= mY -meff,Y(So,S^) -igYo,Wo + gY^fo), (28) 



where wi^.^j^ ^(.Soj.Vq) is the effective mass and Y stands for the A, S 
and H hyperons. For the determination of the gAu, gl.a and gsa 
coupling constants, the following values of the potentials were used 



C/f = -28MeV, C/^"'=-^30MeV, U^'>=- 



IBMeV. 



(29) 



In general the coupling constants gYa can be decomposed into 
two parts gYa=gr7XYr7, whcrc XYa depends on the value of 
potential Uy'\ The coupling of hyperons to the strange meson <j* 
were obtained from the following relations [50] 



(A), 



:2{7 



(A) 



(30) 



The scalar coupling constants are collected in Table 3. 



The symmetry energy 

The TMl parameter set (Table 4) that successfully describes the 
ground state properties of both finite nuclei and infinite nuclear 
matter was supplemented with the mixed nonlinear isoscalar- 
isovector CO — p meson coupling, which provides the additional 
possibility of modifying the high density components of the 
symmetry energy. The remaining nuclear matter ground state 
properties were left unchanged. The symmetry energy is given by 
the relation 



5y^ 



-. + - 



14 



(31) 



+m 



eff,N 



eff,P 



where m^,jj- denotes the effective mass of the p meson in the case 
of non-strange, symmetric nuclear matter. 

The strength of the w — p coupling is set by the Aj^(gQ,gp)^ term 
and requires the adjustment of the gp couphng constant to keep 
the empirical value of the symmetry energy Esym(nb) = 25.68 MeV 
at the baryon density rib, which corresponds to A:ir = 1.15 fm~' 
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Table 3. Scalar meson coupling constants. 



l/^'' (MeV) 


t/i"' (MeV) 


g\„ 




A'E„ 




gZF 


gZi,- 


+30 


-14 


6.169 


3.084 


4.476 


5.626 


1 1 .474 


5.626 


+30 


-18 


6.169 


3.201 


4.476 


5.482 


11.372 


5.482 


+20 


-18 


6.169 


3.201 


4.768 


5.482 


1 1 .372 


5.482 


+10 


-18 


6.169 


3.201 


5.060 


5.482 


11.372 


5.482 


-10 


-18 


6.169 


3.201 


5.644 


5.482 


1 1 .372 


5.482 


-20 


-18 


6.169 


3.201 


5.935 


5.482 


11.372 


5.482 


-30 


-18 


6.169 


3.201 


6.227 


5.482 


1 1 .372 


5.482 


Scalar meson 


coupling constants gYa 


and gY„- [Y 


= A,S,E) calculated for chosen values of the potentia 


s and C/f . 







doi:1 0.1 371 /journal.pone.Ol 06368.t003 



[40]. The parameters Ay together with the adjusted value of the 
parameter gp are presented in Table 5. 

The density dependence of the symmetry energy can be 
expressed by coefiicients that define the slope L and the curvature 
Ksyni of the symmetry energy 



Esyn,(nb)=Esym(no)+L 



3«o 



+ 



2! 



3wo 



(32) 



Equation (32) is a typical low density expansion, higher-order 
terms must be taken into account at suprasaturation densities. The 
performed calculations that were focused mainly on the slope 
parameter of the symmetry energy 



L = 3«( 



5^ £■,,,„,(«/,) 



0- 



(33) 



gave results that agreed with the experimental data (L = 88 + 25 
MeV [51]). These results are collected in Table 5. However, 
current estimation of the symmetry energy parameter based on 
theoretical, experimental and observational results narrows the 
range of L to (40.5-61.9) MeV [52]. 

Results 

The energy density and pressure given by relations (26) and (27) 
define the EoS. Numerical calculations, which were made for the 
TMl parameterisation for both nonstrange and strangeness-rich 
matter, led to the solutions that are shown in Fig. 1. For this 
parameterisation a class of EoSs was obtained. Individual EoS is 
parametrized by the coupling constant Kv> which determines the 
strength of the mixed vector meson interactions. The form of the 
EoSs allows one to compare the differences between various 
models. 

An analysis performed for the nonlinear model with parameter 
A)/, which ranges between 0.014 to 0.017, showed that the 
increase of A[/ produces a stifiFer EoS. For comparison, the EoSs 
for non-strange matter for NL3, TMl and FSUGold parameter- 
isations have been included. Thus, the stiffness of the EoS depends 
on the existence and strength of the mixed vector meson 
interactions, and the nonlinear model described by the Lagrangian 
function (6) makes it possible to construct a much stiffer EoS than 
the one obtained for the TMl -weak model (the abbreviation 
TMl -weak denotes the standard TMl model extended to the fuU 



octet of baryons with two additional meson fields, which were 
introduced in a minimal fashion, to reproduce the hyperon- 
hyperon interaction). 

In the case of nuclear matter an extended isovector sector 
comprises the co — p meson interaction and the parameter Ky sets 
the strength of the 03 — p coupling. This term altered the density 
dependence of the symmetry energy. The standard TMl 
parameterisation without w — p coupling gives as a result very 
stiff form of the symmetry energy. The inclusion of the CD — p 
coupling softens the symmetry energy. The solutions were 
presented in Fig. 2. Calculations were done for rather high values 
of Aj/ = 0.0165 and 0.03. The interaction between o and p 
mesons leads to the solution, which approaches that obtained for 
the AVI 4 and UV14 models with the Urbana VII (UV VII) three 
nucleon potential. For comparison the form of the symmetry 
energy calculated for the UV14 plus TNI model was included 
[53]. 

Interacting baryons are the basic components of the matter of 
neutron stars. The modification of baryon masses that arises from 
baryon interactions with the background nuclear matter is shown 
in Fig. 3. The numerical solutions predicted by equation (18) for 
the fixed value of parameter A(/ = 0.0165 for both the nonlinear 
model and for the TMl -weak model, show reduced effective 
baryon masses. 

The reduction of the nucleon mass in the nonlinear model is 
essentially the same as that obtained in the TMl -weak one. Thus, 
the influence of the nonlinear vector meson couplings on the 
nucleon effective mass is negligible. The effective masses of strange 
baryons in the case of the nonlinear model drop less rapidly than 
the effective masses obtained in the TMl -weak model. The 
behaviour of the baryon effective masses is governed by the density 
dependence of the scalar mean fields, which is presented in Fig. 4. 
The presence of nonlinear couplings between vector mesons 
modifies the density dependence of the strange scalar meson 
leaving the nonstrange scalar meson almost unchanged. 

The in-medium reduction of baryon masses is equivalent to the 
modification of vector meson masses in the meson sector (Fig. 5). 
An analysis of the density dependence of the effective p and (j> 
vector meson masses led to the conclusion that their modification 
was produced by a strong Ay dependence, especially in the high 
density limit. The effective mass of the oj meson is almost 
independent of the value of parameter Ay. 

The stiffness of the EoS is characterised by the incompressibility 
of nuclear matter. In general, incompressibility comprises terms 
resulting from the kinetic pressure of Fermi gas and from the 
potential of the model. Analysing the strange sector of the model 
one can compare the factor that determines the difference between 
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Table 4. The standard TMl parameter set [27]. 





m„ = 511.2 MeV 


= 10.029 


g, = 7.2325 fm - ' 


m„ = 783 MeV 


g„j = 12.614 


gi = 0.6183 


nip = 770 MeV 


& = 9.264 


c, = 71.0375 


doi:l 0.1 371 /journal.pone.Ol 06368.t004 



the strength of the efTective repulsive and attractive forces. The 
scalar meson a* has been introduced in a minimal fashion thus 
f^efi\a' =mnt. The strength of the effective repulsive force between 
strange baryons is mainly altered by the factor 1 / m^ff ip . The 
influence of the Ay parameter on the density dependence of 
l/mi,ff ij, is depicted in Fig. 6. The increase of the parameter Ay 
considerably enhances the strength of the repulsive force in the 
system. 

The properties of asymmetric strangeness-rich matter of 
neutron stars are characterized by the parameters that define 
the strangeness content of the system fs = nsjnb, where ns is the 
strangeness density and the isospin asymmetry fT, = n^jni,, where 
«3 is the isospin density given by the relation 



«3 = ^ hsnB 



(34) 



The density dependence of the asymmetry parameter (22) 
and the strangeness content of the system fs (23) are depicted in 
Fig. 7. The nonlinear model leads to a system with an enhanced 
asymmetry and a considerably reduced strangeness. An increase of 
parameter Ay causes the matter to become more asymmetric. 
Charge neutrality and the condition of p equilibrium 
(p + e^<-^n + Vp) impose constraints on a neutron star composition. 
Assuming that neutrinos are not considered, since their mean free 
path is longer then the star radius, the following relation can be 
obtained: /i^, + ^ip = /i„. 

Additional hadronic states are produced in neutron star 
interiors at sufiicientiy high densities when hyperon in-medium 
energy equals their chemical potential. The higher the density the 
more various hadronic species are expected to populate. In general 
weak reactions for baryons and the corresponding equations for 
chemical potentials can be written in the form 



(35) 



(36) 



where Bi and Bi denote baryons, / and v/ lepton and neutrino of 
the same flavor, whereas and refer to baryon B with baryon 
number bs and charge qs- The above equations create relations 
between chemical potentials of particular hyperons 



(37) 



The presented equilibrium conditions determine all constituents 
of the matter of neutron stars. The concentrations of a particular 
component i = B,l of the matter of a neutron star can be defined as 
Yj = ni/ni,, where «, denotes the density of the component / and 
is the total baryon number density. The density fraction of 
nucleons and leptons as a function of the baryon number density 
for the frxed value of parameter Ay is presented in Fig. 8. The 
important findings concern the concentration of leptons, which are 
more highly populated in the case of the nonlinear model. Fig. 9 
shows how the modification of the vector meson sector alters 
concentrations of strange baryons. The first hyperon that appears 
is A and it is followed by E ^ and 2 ^ . The appearance of 
negatively charged hyperons reduce the concentrations of leptons. 
This stems from the charge neutrality condition. However, the 
initial rapid increase in population of hyperons has been 
suppressed leading to significantly reduced concentration of 
hyperons at sufTicientiy high density. For comparison the results 
obtained for the TMl -weak model have been included. It is 
evident that the additional nonlinear couplings between vector 
mesons modify chemical composition of the neutron star shifting 
the hyperon onset points to higher densities and reduces the 
strangeness content of the system. 

A composition and concentrations of hyperons calculated in the 
nonlinear model can be traced in a chosen configuration of a 
neutron star. The composition of the core of the maximum mass 
configuration is depicted in Fig. 10. The number density of 
hyperons is reduced. However, an interesting feature of this model 
is the abundance of hyperons in the very inner part of the 
neutron star inner core. 

The most essential feature of the model is connected with the 
fact that even in the presence of hyperons the obtained EoS is very 



Table 5. TMl parameter set - isovector sector. 



Av 0 0.014 0.015 0.016 0.0165 0.017 



gp 


9.264 


9.872 


9.937 


10.003 


10.037 


10.071 


L (MeV) 


108.58 


77.52 


75.81 


74.16 


73.36 


72.56 



TMl parameter set with the extended isovector sector. The extension includes the coupling between the co and p mesons: Ay{g^jgp f(cOpp'p-. For each value of the 
parameter Ac the g^ coupling constant has been adjusted to reproduce the empirical value of the symmetry energy Esy„,(nb) = 25. 6S MeV at the baryon density nt 
which corresponds to kf = 1.15//J7"' [40]. 
doi:l 0.1 371/journal.pone.Ol 06368.t005 
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Figure 1. The pressure vs density calculated for selected parameterisations. The stiffest EoSs have been obtained in the case of non- 
strange matter for NL3 and TIVI1 parameter sets. The EoSs calculated for the extended nonlinear TM1 model for different values of the parameter Ky 
form a distinct class with the stiffest EoS obtained for Ky =0.017. The softest EoS for strangeness rich neutron star matter represents the case of the 
standard TM1 model extended by the inclusion of strange mesons which have been introduced in a minimal fashion (TMI-weak model). The EoSs 
calculated for the extended nonlinear model for different values of parameter are located between these two curves. Nonlinear vector meson 
couplings stiffen the EoS which results in a higher value of the maximum mass. 
doi:1 0.1 371/journal.pone.01 06368.g001 



stifT. Global neutron star parameters such as the mass and radius 
and the structure of a neutron star can be determined by the 
equation of hydrostatic equilibrium - the Tolman-Oppenheimer- 
Volkoff (TOV) equation: 



dV{r) _ - G(£{r) + V(r)) {m(r) + Anr^V(r)) 



dr 



2Gm(r) 



250 



200 - 



> 150 



100 - 



50 



A„ = 0 




A„ = 0.0165 




A„ = 0.03 ^/ 




AV14+VII 




UV14+VII 




UV14+TNI 




j^SassS-m — ««:;:: 









0.2 



0.4 



0.6 



0.8 
nb(fm-^) 



1.2 



1.4 



Figure 2. The density dependence of the symmetry energy. The 

density dependence of symmetry energy calculated for TIVI1 parameter 
set [27] for different values of the parameter Ay. For comparison the 
results obtained for the AV14+VII, UV14+VII and UV14+TNI models [53] 
are included. The inclusion of Ay parameter softens the symmetry 
energy and its density dependence resembles that obtained for the 
realistic nuclear models. 
doi:1 0.1 371/journal.pone.01 06368.g002 



dm{r) 
dr 



1400 
1200 
1000 
800 
600 
400 
200 
0 















"""/i;,;. 













































































2 3 4 5 

n, TMI-weak S, TMI-weak ■ 

n, Av=0.0165 Z, Av=0.0165- 

A, TMI-weak H, TMI-weak ■ 
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Figure 3. iVIedium modification of baryon masses. The effective 
baryon masses as a function of relative baryon density calculated for 
the extended nonlinear model, for Aj' = 0.0165 and for the TMI-weak 
model. The extended nonlinear model, in comparison to the TMI-weak 
one, leads to the higher effective strange baryon masses. 
doi:1 0.1 371/journal.pone.01 06368.g003 
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Figure 4. Density dependence of the scalar meson fields. The 

scalar mesons a and a* as a function of baryon density. For comparison 
the results obtained for the TMI-weak model have been included. 
There is a significant difference between the density dependence of the 
strange scalar meson for the TMI-weak and extended nonlinear 
models. The density dependence of scalar meson a is left unchanged 
for both models. 

doi:1 0.1 371 /journal.pone.01 06368.g004 



dn(r) 
dr 



--ATirn(r)\ 1 



2Gm{r) 



-1/2 



(38) 



where m and n denote the enclosed gravitational mass and baryon 
number, respectively, V is the pressure and £ is the total energy 
density. In order to get the numerical solution of equation (38), the 
EoS has to be specified. A correct model of a neutron star is based 
on the assumption that its internal structure is composed of 
separate parts, thus the construction of the mass-radius relation 
requires taking into account additional EoSs that describe the 
matter of the inner and outer crust. For the outer and inner crusts 



Figure 6. The range of hyperon-hyperon interaction. The density 
dependence of the factor l/ni^./j ,j, calculated for the nonlinear model, 
for Ak = 0.014 and 0.0165. The behaviour of this parameter is correlated 
with the density dependence of (j) meson effective mass. 
doi:10.1 371/journal.pone.01 06368.g006 



the EoSs of Baym, Pethick, and Sutherland (BPS) [54] and Baym, 
Bethe and Pethick (BBP) [55] have been used respectively. 

The results obtained for the set of EoSs calculated in this paper 
led to the mass-radius relations and permitted the value of the 
maximum mass to be determined which in a sense can give a 
measure of the impact of particular nonlinear couplings between 
vector mesons. The mass-radius relations obtained for varying 
values of parameter Ay are shown in Fig. 1 1. The higher the value 
of Ay, the higher the maximum mass. 

The influence of the (7^' potentials 

There are stiU significant uncertainties associated with the 
experimental data on the hyperon-nucleus interactions. Thus it is 
reasonable to investigate the effect of the hyperon-nucleus 
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Figure 5. Medium modification of vector meson masses. The 

effective vector meson masses as a function of baryon density 
calculated for the extended nonlinear model. The results for 
Aj/ = 0.014 and 0.0165 have been compared. In the case of oj and p 
mesons the increase of Aj/ parameter results in their higher effective 
masses. While for the effective mass of the ijfi meson there exists a range 
of density where its effective mass decreases. 
doi:1 0.1 371/journal.pone.01 06368.g005 
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Figure 7. The isospin asymmetry and strangeness content of 
the matter of neutron stars. The asymmetry parameter /„ and the 
strangeness content of the system fs as a function of baryon density 
calculated for the extended nonlinear model for different values of the 
parameter Ay. For comparison the result obtained and for the TMI- 
weak model has been included. In the case of the nonlinear extended 
model the matter of the neutron star is less symmetric and with the 
reduced strangeness content. 
doi:1 0.1 371/journal.pone.01 06368.g007 
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u = ni, / no 



Figure 8. Nucleon and lepton concentrations in tKie matter of 
neutron stars. Relative concentrations of nucleons and leptons 
calculated for the extended nonlinear model, for the fixed value of 
the parameter Af^ = 0.0165. For comparison the results obtained for the 
TMI-weak model have been included. In the case of the TMl-weak 
model the lepton concentrations are reduced whereas the differences 
in neutron and proton concentrations are rather small. 
doi:1 0.1 371 /journal.pone.Ol 06368.g008 

potential U^y^ on the obtained results [56]. Particular attention 
was paid to tlie dependence of the EoS on the S -nucleus potential 
Uf. 

Detailed calculations were done for tlie selected values of the 
potential, assuming its both attractive and repulsive 
character. Calculations performed for the extended nonlinear 
model were resulted in a sequence of EoSs (Fig. 12). The stififest 
one was obtained for the repulsive potential = 30 MeV). For 
comparison the EoS obtained for TMl-weak model with the 
repulsive U'^^ = 30 MeV potential was included. In order to study 
the influence of the U^^'' potential a model with an exaggerated 
value of the potential L'^'^' = — 1 MeV was examined. The 




Figure 9. Concentrations of strange baryons in the matter of 
neutron stars. Relative concentrations of strange baryons calculated 
for the extended nonlinear model, for the fixed value of parameter 
A|/ = 0.0165. For comparison the results obtained for the TMl-weak 
model have been included. The extended nonlinear model leads to the 
reduced number density of A and hyperons. Characteristic feature 
of this model is the appearance of hyperons. 
doi:1 0.1 371/journal.pone.01 06368.g009 




r( km ) 

Figure 10. Chemical composition of the maximum mass 
configuration. The particle fraction y, =«,/«;, as a function of the 
radius of the neutron star, for the maximum mass configuration. 
Calculations have been done on the basis of the extended nonlinear 
model with Ap/ = 0.0165. The dotted vertical lines point the threshold 
densities for particular hyperons. 
doi:10.1371/journal.pone.0106368.g010 

resulting change in the EoS is negligible. A similar conclusion can 
be drawn by examining the changes caused by the reduction in the 
potential C/l^. Calculations were done for two chosen values of 
the t/i^' potential: -14 MeV and -18 MeV. 

Taking into account the value of the maximum mass obtained 
for a given EoS the presented results support the conclusion that 
the value of the potential is not a factor that decisively influences 
the form of the EoS. A similar conclusion can be obtained by 
analysing the mass-radius relation (Fig. 13). The most promising 
results were obtained for the repulsive = 30 MeV potential 
which within the considered model leads to the highest value of the 
maximum mass. 

Hyperons are distributed in the very inner part of a neutron star 
core therefore a schematic cross-section shows a hyperon core 
inside of a strangeness-rich neutron star. Information concerning 




R ( km ) 



Figure 11. The mass-radius relations. The mass-radius relations 
calculated for the nonlinear model for different values of parameter Ay- 
The results have been obtained for both the TM1 parameterisation in 
the case when the matter of the neutron star comprises only nucleons 
and for the FSUGold parameter set. 
doi:1 0.1 371/journal.pone.01 06368.g01 1 
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Figure 12. The influence of the [4 potential on the EoS. The 

pressure vs density calculated for different values of the t/^'^' potential 
(attractive and repulsive), the case of a very weak A — A interaction 
(L'J^'^' = 1 IVIeV) has been also included. For comparison the EoS 
calculated for the TMI-weak model has been shown. In the case of the 
extended nonlinear model the higher value of the repulsive 
potential leads to the stiffer EoS. 
doi:1 0.1 371 /journal.pone.01 06368.g01 2 



properties of the internal structure of a neutron star are 
summarised in Fig. 14 and 15. These figures depict the 
dependence of both neutron stars and their hyperon cores on 

the C/j'^' potential and strictly speaking on the value of x-za 
parameter and offer indications for interpreting the results of 
numerical calculations. Calculations were done for maximum 
mass configurations. The mass of the star increases when U'^^ 
potential becomes more repulsive whereas the mass of the hyperon 
core decreases reaching a minimum value for the parameter 
•XE(t~0.51 (Fig. 14). Similar behaviour has the radius of the 
maximum mass configuration. Results are presented in Fig. 15. 
One of the most important difference between the model 



Figure 14. The influence of the IJ\_ potential on the internal 
structure of the maximum mass configuration. The dotted line 
shows the dependence of the value of the maximum mass Ms 
obtained for the given EoS on the parameter ax„. The parameter xy_„ 
carries information about the hyperon-nucleon interaction. Solid line 
depicts similar relations for the mass of the hyperon core . 
doi:1 0.1 371/journal.pone.01 06368.g01 4 

the neutron star matter. The hyperon onset points depend on the 
character of the t/j^' potential. Calculations performed on the 
basis of the extended nonlinear model indicate that the 
enhancement of the attractive C/^'^' potential shifts the hyperon 
onset points to lower densities, whereas the change of the repulsive 
potential leaves the hyperon onset points almost unchanged 
(Fig. 16). In the case of attractive potential the first hyperons that 
appear in neutron star matter are 2^ hyperons and they are 
followed by A hyperons. Such a scheme of the emergence of 
hyperons is changed in the nonlinear model with the repulsive 
t/^''' = 30 MeV potential. The first hyperons that appear are A 
hyperons and successively and hyperons. The relative 
concentrations of hyperons calculated for both attractive and 

m 



constructed for the attractive C/f > = - 30 MeV and the repulsive repulsive potential are presented in Fig. 1 7. The results of 



f/j^' = 30 MeV potential is related to the chemical composition of 



numerical calculations obtained for different values of the potential 
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Figure 13. The influence of the {7^ potential on the mass- 
radius relation. The mass - radius relations calculated for EoSs 
presented in Fig. 12. In the case of the extended nonlinear model the 
higher value of the repulsive C^"^' potential leads to the higher value of 
the neutron star maximum mass. 
doi:1 0.1 371/journal.pone.01 06368.g01 3 
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Figure 15. Neutron star and their hyperon core radii as 
functions of the hyperon-nucleon interaction. The dotted line 
shows the dependence of the radius of the maximum mass 
configurations Rs obtained for the given EoS on the value of parameter 
x^„. The solid line depicts similar relations for the radii of the hyperon 
core Rk^. 

doi:1 0.1 371/journal.pone.01 06368.g01 5 
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Figure 16. Hyperon onset points for the maximum mass 
configurations. The dependence of the hyperon onset points on 
the value of the parameter x^„. 
doi:1 0.1 371 /journal.pone.01 06368.g01 6 

J7^'^' were compiled in Table 6. These results for a given EoS 
include the value of the maximum mass, the radius of the 
maximum mass configuration, the radius and mass of the hyperon 
core and the onset point of hyperons. For comparison the model 
with the U^^^ = - 1 MeV was included. 

The problem connected with the impact of the [/'j^' and U\P 
potentials on the gross parameters of a hyperon-rich neutron star 
has been developed so far. Solution of the TOV equations (38) can 
also provide information on the influence of the hyperon-nucleon 
and hyperon-hyperon interactions on the internal structure of a 
given star. The purpose of such an analysis is to find similarities 
and cUfferences in the internal structures of neutron stars whose 
models were obtained for significantly dilferent equations of state. 
Comparison of the internal structure of stars was carried out for 
the configuration of the maximum mass achievable for the given 
model, for the cases with different values of the U^'' potential. 




0 1 2 3 4 5 6 7 

u = nt,/no 



Figure 17. Hyperon concentrations for opposite signs of tKie 
potential. Comparison of the relative concentrations of hyperons 
calculated for the extended nonlinear model for the repulsive and 
attractive V^'^ potential. In the case of the repulsive potential the first 
hyperon that appears is A and it is followed by and E^, for the 
attractive potential hyperons appear in the following order: E^, A and 

doi:1 0.1 371 /journal.pone.01 06368.g01 7 



Figure 18. The profiles of the asymmetry of the matter of 
neutron stars. The asymmetry of the matter of neutron stars for the 
maximum mass configurations. Calculations have been done for the 
extended nonlinear model for the fixed value of the parameter 
Ak = 0.0165. Profiles of the asymmetry parameter /„ have been 
estimated for the U^' potential ranging from +30 iVleV to —30 
IVleV. The asymmetry derived for the TlVll-weak model for the repulsive 

=30 iVleV potential is also included. 
doi:1 0.1 371/journal.pone.01 06368.g01 8 

The results that relate to the change of asymmetry parameter /„ 
and the strangeness content of the star fg are presented in Fig. 1 8 
and 19 respectively. In both cases the obtained results show the 
existence of significant differences between the TM 1 -weak model 
and the extended nonlinear model. Considering the asymmetry 
parameter fa the largest differences occur in the outer part of the 
hyperon core where TMI-weak model leads to the enhanced 
asymmetry of the system. Comparing the profiles of the fs 
parameter it is evident that the very inner part of the hyperon core 
is characterized by markedly enhanced strangeness content of the 
matter modelled on the basis of TMI-weak parameter set in 
contradiction to the results obtained for the outer part of the 
hyperon core. 

The differences in the radial dependence of the parameters fa 
and fs are connected with the chemical composition of the matter 
of neutron stars especially with the reduced abundance of 
hyperons in the extended nonlinear model. The influence of 
U^'' potential on the relative concentrations of A hyperons for the 
maximum mass configuration is depicted in Fig. 20. A hyperons 
are the most abundant in the case of repulsive potential = 30 
MeV. For lower value of the repulsive potential concentrations of 
these hyperons are reduced. The lowest concentrations of A 
hyperons occur for the attractive potential = — 30 MeV. 

Calculations obtained for the TMI-weak model for the repulsive 
' = 30 MeV potential give as a result much higher concentra- 
tion of A hyperon. In the case of S hyperons the change of the 
potential produces the opposite effect to that which is 
obtained for A hyperon. The highest abundance of S ^ hyperons 
(Fig. 2 1) is obtained for the attractive potential U^^ = — 30 MeV. 
For other cases of the t^f ' potential the concentration of 
hyperons is dropping reaching the lowest value for the repulsive 
potential U^^'* = 30 MeV. In the case of the TMI-weak model 
with the repulsive U^'' = 30 MeV potential the relative concen- 
tration of hyperons vanishes. hyperons (Fig. 22) occur 
inside the hyperon core for suflficiendy large values of the attractive 
potential (-20 MeV and -30 MeV). hyperons (Fig. 22) 
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Table 6. The influence of the Uy potential on neutron star parameters. 





£4'^' (MeV) 


f/f (MaV) 


M. (Mq) 


R, (km) 




Rlu (km) 


nb.aia,i {ft" ^) 


+30 


-1 


2.036 


10.47 


1.620 


8.32 


0.348 


+30 


-5 


2.032 


10.45 


1.617 


8.31 


0.348 


+20 


-5 


2.027 


10.44 


1.613 


8.29 


0.348 


+10 


-5 


2.020 


10.31 


1.609 


8.27 


0.348 


-10 


-5 


2.008 


10.23 


1.679 


8.38 


0.325 


-20 


-5 


2.003 


10.18 


1.720 


8.51 


0.305 


-30 


-5 


2.001 


10.04 


1.768 


8.59 


0.288 



Neutron star maximum masses and corresponding radii calculated for the chosen values of the potentials and ^7J^ . A/fa. and J?/„. denote the mass and radius of a 
hyperon core for a given maximum mass configuration, n/,_ai,si-: denotes the baryon density at which hyperons appear in the system (hyperon onset point). 
doi:l 0.1 371 /journal.pone.Ol 06368.t006 



appear in the hyperon core only in the case of the repulsive 
potential. The higher the value of the repulsive potential, the 
higher the abundance of hyperons. 

Discussion 

Recent observations of the binary millisecond pulsars give in the 
result the value of neutron star masses of the order of two solar 
masses. Theoretical models that describe strangeness-rich matter 
of a neutron star lead to much softer EoSs than those constructed 
for a neutron star composed of only nucleons. This significant 
softening of the EoS results in the low value of the maximum 
neutron star mass achievable in a theoretical model of a neutron 
star that involve hyperons and in general, leads to the 
inconsistency between theoretical models and observations. 

One of the possible solution to this problem, is the construction 
of a theoretical model of the strangeness-rich matter of a neutron 
star that will result in the emergence of extra repulsive force in the 
strange sector of the system. This in turn stiffens the EoS and leads 
to the higher value of the maximum mass of a neutron star. 



Q g ; TM1-weak 

U'^j = + 30 MeV 

0.5 ■ '\ u'^j = + 20MeV 

U'^j = + 1 0 MeV 
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Figure 19. The strangeness content of the matter of the 
maximum mass configuration. Results obtained on the basis of the 
extended nonlinear model for the fixed value of the parameter 
Ak = 0.0165. Profiles of the parameter fs have been estimated for the 
potential ranging from +30 MeV to —30 MeV. The strangeness 
content derived for the TMI-weak model for the repulsive f^'^' = 30 
MeV potential is markedly different then those for the nonlinear model. 
In the case of TMI-weak model the hyperon core is more compact but 
with substantially enhanced hyperon content. 
doi:1 0.1 371 /journal.pone.Ol 06368.g01 9 



The model that was used in this paper for the description of the 
matter of a neutron star is based on the nonlinear realization of the 
chiral symmetry. The results were obtained under the assumption 
that nucleons do not couple to meson (j). Alternative analysis of the 
behavior of the EoS based on 5'[/(3) symmetry for different 
relativistic models have been already done and the results are 
included in [57-59]. Whereas the analysis of the EoS of 
hypernuclear matter within a relativistic density functional theory 
with density-dependent couplings has been included in [60]. 

An essential characteristic of the model that was used in this 
paper is the diverse vector meson sector which comprises different 
vector meson couphngs. The performed analysis concerns the issue 
of the influence of the nonlinear co, p and ^ vector meson 
interactions on the properties of asymmetric strangeness-rich 
nuclear matter and consequendy on neutron star parameters. 

As a result, a special class of EoSs was obtained. Particular EoSs 
are characterised by the value of parameter Ky, which determines 
the strength of the vector meson couplings. This very special form 
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Figure 20. The influence of {7^^' potential on the A hyperon 
fractions for the maximum mass configuration. The figure 
depicts the change in the concentration of A hyperons caused by 
different values of the U^' potential. A hyperons are the most 
abundant in the case of repulsive potential U^''^30 MeV. For lower 
value of the repulsive potential concentrations of these hyperons are 
reduced. The lowest concentrations of A hyperons occur for the 
attractive potential (7^'^' = —30 MeV. Calculations obtained for the 
TMI-weak model for the repulsive U^'' = 30 MeV potential give as a 
result much higher concentration of A hyperon. 
doi:1 0.1 371/journal.pone.01 06368.g020 
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U"j. = + 30 MeV 
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Figure 21. The influence of potential on the £ hyperon 
fractions for the maximum mass configuration. In the case of E 
hyperons the change of the potential produces the opposite 

effect to that which is obtained for A hyperon. The highest abundance 
of S~ hyperons is obtained for the attractive potential U^'=—30 
MeV. For other cases of the potential the concentration of 
hyperons is dropping reaching the lowest value for the repulsive 
potential t/^"' = 30 MeV. In the case of the TMI-weak model with the 
repulsive U^^ = 30 MeV potential the relative concentration of 
hyperons vanishes. 

doi:10.1371/journal.pone.0106368.g021 

of the considered model of matter of a neutron star with hyperons 
gives an EoS that is much stiffer than that obtained using the 
standard TMI-weak model in which hyperons are introduced in a 
minimal fashion. 

Obtained numerical solutions allow one to analyse the impact of 
the nonlinear vector meson couplings on the density dependence 
of the scalar and vector meson fields. It was shown that under the 
condition of the presence of nonlinear vector meson couplings the 
in-medium properties of baryons and mesons were changed. 
Especially there is a considerable modification of the effective 
baryon and vector meson masses. 

The nonlinear vector meson couplings modify both the 
asymmetry and strangeness content of the system and therefore 
lead to a model with a reduced strangeness and an enhanced 
asymmetry. The stiffness of the EoS is directly related to the in- 
medium properties of the matter of a neutron star. It depends on 
value of the effective baryon and meson masses, which modify the 
compressibility of the matter of a neutron star and the range of 
interactions especially in the strange sector. The results of the 
analysis performed in the framework of the nonlinear model have 
shown that the properties of neutron stars are significantly altered 
by the presence of hyperons. 

The analysis of the dependence of neutron star parameters on 
the strength of hyperon-nucleon interaction has been already done 
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linear models. The reduction of the hyperon population in the 
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is presented. In this particular model the core of a neutron star 
reveal a hyperon inner core with the reduced A and hyperon 
concentrations but with the enhanced population of 2 ^ hyperons. 
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